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ABSTRACT 

We present the X-ray Luminosity Function (XLF) of the ROSAT Deep Cluster 
Survey (RDCS) sample over the redshift range 0.05-0.8. Our results are derived from 
a complete flux- limited subsample of 70 galaxy clusters, representing the brightest half 
of the total sample, which have been spectroscopically identified down to the flux limit 
of 4 x 10 -14 erg cm~ 2 s _1 (0.5-2.0 keV) and have been selected via a serendipitous 
search in ROSAT-PSPC pointed observations. The redshift baseline is large enough 
that evolutionary effects can be studied within the sample. The local XLF {z < 0.25) 
is found to be in excellent agreement with previous determinations using the ROSAT 
All-Sky Survey data. The XLF at higher redshifts, when combined with the deepest 
number counts constructed to date (/ > 2 x 10 -14 erg cm -2 s _1 ), reveal no significant 
evolution at least out to z = 0.8, over a luminosity range 2 x 10 42 - 3 x 10 44 erg s 1 in 
the [0.5-2 keV] band. These findings extend the study of cluster evolution to the highest 
redshifts and the faintest fluxes probed so far in X-ray surveys. They complement and 
do not necessarily conflict with those of the Einstein Extended Medium Sensitivity 
Survey, leaving the possibility of negative evolution of the brightest end of the XLF at 
high redshifts. 

Subject headings: galaxies: clusters: general — X-rays: general — cosmology: 
observations 
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1. Introduction 

Large homogeneous samples of galaxy clusters, spanning a wide range of redshifts, have long 
been considered powerful tools to study the evolution of the large scale structures in the Universe. 
Unfortunately, finding clusters at cosmologically interesting lookback times, i.e. z > 0.5, not to 
mention defining a complete sample, is a time consuming and difficult task. As a result, the 
much-needed observational constraints for theories of cluster formation have not been forthcoming. 

In an effort to remedy this situation, we have embarked on a project, the ROSAT Deep Cluster 
Survey (RDCS), aimed at constructing a large homogeneous sample of distant galaxy clusters 
selected solely on the basis of their X-ray properties from ROSAT-PSPC pointed observations. 



Initial results have been presented in Rosati et al. 1995 (R95). Pure X-ray selection leads to a 



selection function which can be modelled in a relatively straightforward way, being essentially that 
of a flux limited sample. As a result, distribution functions of observables, such as the number 
counts, the X-ray Luminosity Function (XLF) and the redshift distribution, can be directly 
compared with theories of structure formation. Until very recently, the only available sample 
of X-ray selected clusters at high redshifts was the one compiled from the Einstein Extended 
Medium Sensitivity Survey (EMSS) flGioia et al. 1990| , |Henry et al. 1992] (H92)). From a complete 



subsample of 67 clusters at 0.14 < z < 0.6, the authors found evidence of a statistically significant 
negative evolution of the XLF at z ^ 0.3 and ^x[o.3— 3.5] keV ^ 5 x 10 erg s _1 . Earlier claims of 
strong evolution at lower redshifts ( Edge et al. 1990| ) have been ruled out by a large compilation of 



X-ray clusters from the ROSAT All-Sky Survey (RASS) (|Ebeling et al. 1997| ). More recently, the 



EMSS findings have been challenged by a re-analysis of the same sample supplemented by ROSAT 
data (Nichol et al. 1997). Covering a very large solid angle (735 deg 2 ), but being relatively shallow, 



the EMSS sample can probe mostly the bright end of the XLF and can provide weak constraints 
on the evolution of the XLF at higher redshifts since only 6 EMSS clusters lie at z > 0.5. On the 
other hand, the RDCS reaches considerably fainter fluxes, but covers a much smaller solid angle 
and therefore probes a complementary region in the Lx — z plane to the EMSS. Surveys similar to 
the RDCS, utilizing the PSPC archival data and somewhat different selection techniques, are now 
well underway (RIXOS QCastander et al. 1995|) , SHARC (|Collins et al. 1997| , purke et al.^99p , 



WARPS ( gcharf et al.~1997b ). 

In this Letter, we present the XLF constructed from a flux-limited, spectroscopically 
confirmed subsample of the RDCS. This subsample is already as large as the EMSS sample and 
covers a wider redshift range, from the local universe (z=0.05) out to z ~ 0.8. Therefore, we can 
compare our low-redshift luminosity function with recent determinations from cluster samples 
compiled from the ROSAT All-Sky Survey (Burns et al. 199£, Ebeling et al. 19"97| ). Furthermore, 



we can compare, for the first time, the space density of high redshift clusters with the findings of 
the EMSS and extend the study of cluster evolution to even higher redshifts. We adopt Hq = 50 
km s _1 Mpc -1 and qo = 0.5 througout. 
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2. The Cluster Sample 

A full discussion of the analysis of the X-ray data and the selection technique for the RDCS 
sample is presented in R95. Cluster candidates are selected from a serendipitous search for 
extended X-ray sources in deep pointed observations drawn from the ROSAT-PSPC archive. A 
wavelet-based technique, which is particularly efficient at reducing the effects of confusion and is 
not biased against low surface brightness features, is used to detect and characterize X-ray sources. 
A control sample of about 5000 point sources is used to establish the mean and the variance of 
the PSF at different off-axis angles (8), thereby assessing the statistical significance for a source 
to be extended. Cluster candidates are selected as those sources extended at 99% confidence level 
and having a flux /_i 4 = /x[0.5 — 2.0]keV/(10~ 14 erg cm~ 2 s _1 ) > 1 and 8 < 18 arcmin. This 
selection technique yielded 160 cluster candidates over an area of about 48 deg 2 , drawn from 170 
X-ray fields scattered across the two galactic caps (\b\ > 20°). 

The completeness flux limit of the survey is determined by the flux level at which extended 
and point-like emission can be reliably distinguished. In addition to the source flux, this critically 
depends on the off-axis angle 8 within which the candidates are selected, due to PSF degradation. 
From a consideration of surface brightness dimming alone, one expects the fraction of clusters 
which are unresolved by the PSPC at high redshifts to increase at faint flux levels and at high 
off- axis angles. These effects are analyzed in detail elsewhere (Rosati & Delia Ceca, in preparation). 
The selection function is modelled and quantified using a combination of simulations, control 
samples of known distant clusters and by studying cluster number counts as a function of (fx, 8). 
This analysis allows two statistically complete and independent samples to be defined: one deeper 
(Sample A), covering ~33 deg 2 (/-14 > 1,8 < 15') and one shallower (Sample B) with /_i4 > 6 
and 15' < 8 < 18'. By combining these distinct samples the surveyed area increases to ~48 deg 2 . 
We have chosen a higher flux limit for Sample B since it becomes significantly incomplete for 
/-14 < 5. The incompleteness of Sample A is measured to be negligible for f-14 > 4, is about 10% 
at /_i4 ~ 2 and drops to ~ 20-25% at 1 < /_i 4 < 2. 

The sky coverage £1 of the survey, as emphasized in R95, also depends on the intrinsic angular 
size (0) of the source for a given flux, i.e. £1 = f2(/, ©). Here, for the sake of simplicity, we have 
reduced this bivariate function to a function of the flux only, by choosing a value equal to the 
median angular size of the clusters in the whole sample after deconvolving the PSF (0 ~ 58"). 
We have verified that the use of this effective sky coverage is equivalent to the general method 
described in R95. The function Q = Cl(f, 0) for Sample A and its extension (A + B) is shown in 
fig.l. 

3. Optical Follow-up Observations 

In order to identify these cluster candidates, we have undertaken a large optical follow-up 
program, consisting of deep imaging using the KPNO 4m and 2.1m and the CTIO 4m and 
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1.5m, and multislit spectroscopy carried out with KPNO 4m and the ESO 3.6m. The imaging 
survey in V and I bands is now 90% complete and has revealed a very high success rate of 
identification, with about 115 new clusters (or groups) confirmed to date, i.e. displaying a 
significant overdensity of galaxies around the peak of the X-ray emission. More importantly, the 
spectroscopic follow-up work has secured 95 cluster redshifts so far, spanning the range 0.045-0.83. 
Most of the spectroscopic identifications are based on more than 3 redshifts within 1500 km/s, 
with a median value of 5 members per cluster. A significant fraction of the newly discovered 
clusters lie at high redshift, about one-third at z > 0.4 and a quarter at z > 0.5, making the 
RDCS the largest sample of spectroscopically confirmed distant clusters compiled to date. Details 
on the optical observations of the whole RDCS catalogue will be published elsewhere. 



4. The X-ray Luminosity Function 

The imaging and spectroscopic identification of the cluster candidates in Sample A is 97% 
complete down to /_i4 = 4, hence a complete flux limited subsample of 70 clusters with measured 
redshifts can be defined^. We split the sample in three redshift shells, [0.045-0.25], [0.25-0.5], 
[0.5-0.85], where we have respectively 33, 23, 14 clusters. The luminosities in the observed [0.5-2.0] 
keV band range from 3.5 x 10 44 to 1 x 10 42 erg s _1 and thus probe the region of moderately rich 
clusters to poor groups. A non-parametric representation of the XLF has been obtained using 
the \/V a method of Avni and Bahcall (1980). The application of this procedure for deriving the 
XLF in different redshift shells is fully described in Maccacaro et al., 1991. For each luminosity 
bin containing n objects, the differential XLF is computed as, dN(L)/dL = Ya=i l/X^Si ^L), 
where Vs t represents the total search volume for object i and AL is the width of the luminosity 
bin. Bins of equal logarithmic width ALog L = 0.3 have been used. The corresponding 68% error 
bars have been determined using Poissonian statistics (following Wolter et al. 1994). A power law 
spectrum with energy index 0.5, which well approximates a Raymond-Smith spectrum over a large 
temperature range in the [0.5-2.0] keV band, has been used for computing k-corrections. Count 
rates have been converted to fluxes according to R95. 

The XLF derived is shown in fig. 2. We point out the excellent agreement between the local 
XLF of the RDCS, i.e. in the lowest redshift bin, and two independent determinations of the local 
XLF from [Ebeling et al. 1997| and Burns et al. 1996| (B96). Both these surveys use RASS data 



but completely different selection techniques. The Bright Cluster Sample (BCS) of Ebeling et al 



1997| is an X-ray flux limited sample out to z = 0.3 thus covering the intermediate and bright end 
of the XLF, whereas the B96 sample is an optically selected, volume complete, sample of nearby 
groups and poor clusters {z < 0.15) which probes mostly the faint end, down to luminosities well 



6 To date, only two candidates have not been identified down to the flux limit of f-u = 4. These are extended X- 
ray sources which may be either a blend of multiple sources, for which we did not find an obvious optical counterpart, 
or high redshift clusters, for which deeper imaging is required. 
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below 10 42 erg s -1 . The RDCS XLF in the first redshift shell is well fit by a power law in the 
form n{La) = KL^, where L44 is the X-ray luminosity in the the rest frame band [0.5-2.0]keV 
in units of 10 44 erg s _1 and K is in units of 10 -7 Mpc — 3 [L44]( a— •'J. A maximum-likelihood fit 
yields: a = 1.83 ± 0.15 and K = 3.26 ± 0.57. These parameters compare favorably with those of 
B96, a = 1.71 ± 0.19 and K = 3.50 ± 0.77 (after rescaling to H = 50 km s" 1 Mpc" 1 ), and the 
power law part of the Schechter XLF of the BCS sample, a = 1.85 ± 0.09 and K = 3.32(±g;||). A 
different cluster sample, analogous to the BCS, compiled from the RASS in southern sky provides 



similar results ( De Grandi, 1996| ). The excellent agreement found between these independent 



determinations is by no means trivial since the RDCS is the only sample of the four whose selection 
is not driven by any optical information, but purely by the X-ray properties of clusters. Indeed, it 
suggests that systematic effects in local cluster samples seem to be now well under control. 

The inspection of the RDCS XLF in the higher redshift bins does not show any significant 
evolution, at least out to z ~ 0.8, over the probed luminosity range. Based on the number of 
high redshift clusters in the SHARC-South sample |Collins et al. 1997 also found no evidence of 



significant negative evolution. More recently, Burke et al. 1997 have contructed the XLF from the 



same small SHARC-S sample, 16 clusters in the range 0.3 < z < 0.7, and have reached similar 
conclusions to the ones reported in this work. The large fraction of high redshift clusters in the 
RDCS cannot be reconciled with the results of the RIXOS survey QCastander et al. 1995|) . A factor 
of 4 more RDCS clusters have been identified at z > 0.4 compared to the RIXOS survey down to 
the same flux limit (/-14 = 3) and over the same solid angle. This difference may partly be due to 
incompleteness or small number statistics in RIXOS, but also may reflect an overestimate of the 
sky coverage of the RIXOS survey at low fluxes O) is assumed constant) which could mimic 
a dearth of clusters at high redshifts ( Rosati fc Delia Ceca, 199li| ). 



Given its relatively small surveyed area, the RDCS cannot probe the XLF above 
Lx — 3 x 10 44 erg s . This is the luminosity range where the EMSS sample finds evidence of 
negative evolution, i.e. a steepening of the high end of the XLF at z ^ 0.3 ( |Gioia et al. 1990| ). A 
direct comparison of the cluster volume densities of the RDCS and EMSS samples is illustrated 
in fig. 3, where we recompute the XLF of the RDCS sample in the redshift shell [0.3-0.6] (24 
objects), the highest redshift shell of the EMSS sample. The XLF data points of the EMSS 
sample (open and filled squares) have been converted to the 0.5-2.0 keV band by recomputing 
the XLF using the sample, the sky coverage and the method of H92, updating a few redshifts 



and identifications according to Gioia & Luppino 1994. A power law spectral model with energy 



index 0.5 has been used to convert fluxes from 0.3-3.5 to 0.5-2.0 keV band. We have verified that 
the same XLF data points of H92 are obtained within 10% in the original 0.3-3.5 keV band. We 
note the excellent agreement between the RDCS and the EMSS samples in the lowest redshift bin, 
where the respective XLFs overlap. More importantly, the EMSS XLF in the high redshift shell is 
consistent within the errors with the RDCS XLF over the luminosity range 8 x 10 43 - 3 x 10 44 . 
The RDCS however, cannot follow the drop of the EMSS at the bright end of the XLF. 
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5. Cluster Number Counts 

In addition to the spectroscopically confirmed sample, a larger number of cluster candidates 
have been identified from imaging data. We can now compute the cluster number counts from a 
sample of 130 candidates drawn from Sample A+B down /_i4 = 2, 90% of which have positive 
identification. We do not attempt here to extend this analysis below /_i4 = 2 since the optical 
identification is still largely incomplete at such low fluxes. The resulting cluster log iV-log S, which 
significantly improves the one of R95, is shown in fig. 4 along with previous determinations. The 
error bars take into account the statistical uncertainties as discussed in R95 and the systematics 
due to flux measurements. The latter typically increase the upper error bar of N(>S) by +15% 
at high fluxes as a result of the possible lost flux in the evaluation of the wavelet algorithm. The 
faintest three data points have been corrected for the small 10% incompleteness. These surface 
densities are in good agreement with those of the WARPS survey (Jones et al., 1997) above their 
flux limit /_i4 = 6. A maximum likelihood fit to the differential number counts with a double 
power law model, for s > S B , KiS% 2 ~ ai] 1 S~ a2 for s < S B , yields: «i = 2.22 ± 0.15, 

a 2 = 1.91 ± 0.18, S B = 1.03 ± 0.05, Ki = 1.10 ± 0.15 deg" 2 , where fluxes are in units of 
10 -13 erg cm -2 s -1 . This curve can be integrated to provide the cluster contribution to the soft 
X-ray background, Icxb{0-5 — 2keV) ~ 1.5 keV cm _2 s _1 sr _1 down to /-14 = 1, and is found to 
account for about 10% of the observed soft extragalactic X-ray background. 

Given the accurate determination of the local XLF (fig. 2), which we parametrize as the BCS 
Schechter XLF with the RDCS faint slope, we can now compare our cluster logA^-logS with 
the predicted number counts for a non-evolving XLF. In fig. 4 we show that these predictions, 
obtained by integrating the local XLF over the probed redshift and luminosity range for two 
values of qo, are very close to the observed counts. We note that these number counts are 
drawn from a sample which is a factor 2 deeper than the one used to compute the XLF. The 
logA^-logS at low fluxes is very sensitive to the faint end slope of the XLF, as well as its 
evolution, which is not directly probed by the RDCS XLF in the two highest redshift bins. A 
significant steepening (or flattening) (|Aa| > 0.2) of the faint end (Lx < 10 43 erg s _1 ) of the XLF 
at z > 0.3 would lead to an overprediction (or underprediction) of the observed faint number 
counts outside the error bars, for any value of < qo < 0.5. This further strengthens the evidence 
that the XLF does not evolve, within the present uncertainties, over a wide range of luminosities 
(2 x 10 42 <^ Lx{ erg s~ 4 ) ^ 3 x 10 44 ). The implications that these faint number counts have for 



CDM models of cluster formation have been explored by Kitayama Sz Suto (1997) and Mathiesen 
fc Evrard (1997). 



6. Conclusions 



The large size of the RDCS sample and its wide redshift baseline allow us to provide an 
accurate assessment of the evolution of the space density of galaxy clusters over a considerable 
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range in X-ray luminosities and out to the highest redshifts probed so far in X-ray surveys. We 
find an excellent agreement with independent determinations of the local XLF. By combining the 
RDCS XLF with an improved determination of faint cluster number counts, we find that the XLF 
does not evolve significantly at z <J 0.8 in the luminosity range (2 x 10 42 — 3 x 10 erg s _1 ). Our 
results extend and complement those of the EMSS, further restricting the possibility of negative 
evolution in the X-ray cluster population at z ^ 0.3 to only those clusters with X-ray luminosities 
in excess of the local L* Q 5 _ 2 ]k e v — ^ x x ^ 44 er § s_1 - 
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S [0.5-2.0 kev] erg cm 2 s 1 

Fig. 1. — The effective sky coverage of the RDCS for Sample A (solid line) and its extension, 
Sample A + Sample B (dashed line). The thick solid line denotes the spectroscopically confirmed 
subsample used for computing the XLF. 
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L x (0.5-2.0 keV) erg s" 1 



Fig. 2. — The X-ray cluster luminosity function for the RDCS sample in three redshift shells. 
Independent determinations of the local XLF from ROSAT All-Sky Survey data are also shown 
(dashed and solid lines). 
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L x (0.5-2.0 keV) erg s" 1 



Fig. 3. — Comparison of the luminosty functions for the RDCS and the EMSS samples at low and 
high redshifts. 
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Fig. 4. — Observed cluster cumulative number counts for the RDCS sample (Sample A + B) 
and previous determinations. The no evolution curves were computed by integrating the best 
determination of the local XLF (see text) over the luminosity range 1 x 10 42 — 5 x 10 45 erg s _1 out 
to z = 1.1. 



